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Abstract. 

We study the expansion of a Bose-Einstein condensate trapped in a combined 

' optical-lattice and axially-symmetric harmonic potentials using the numerical solution 

^ ■ of the mean-field Gross-Pitaevskii equation. First, we consider the expansion of such 

1 > , 

, a condensate under the action of the optical-lattice potential alone. In this case the 

result of numerical simulation for the axial and radial sizes during expansion is in 
agreement with two experiments by Morsch et al. [2002 Phys. Rev. A 66 021601 (R) 
and 2003 Laser Phys. 13 594]. Finally, we consider the expansion under the action of 
Q ■ the harmonic potential alone. In this case the oscillation and the disappearance and 

O | , revival of the resultant interference pattern is in agreement with the experiment by 

Miillcr et al. [2003 J. Opt. B: Quantum Semiclass. Opt. 5 538]. 
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O ■ 1- Introduction 

OO 

^ ■ The experimental loading of a cigar-shaped Bose-Einstein condensate (BEC) in both 
one- jHI2] and three-dimensional [3] optical lattice potentials generated by standing- 
wave laser fields has initiated a new class of investigations. The optical-lattice potential 
has been used in the study of matter- wave interference pQ , of oscillating atomic current 
in a one dimensional array of Josephson junctions [HE], of Bloch oscillation and 
Landau Zenner tunneling [S], of superfluid- insulator classical [7] and quantum t 3.j phase 
O ■ transitions, and in the generation of matter- wave bright soliton [8 . 

As the size of the condensates in above experiments is much too small, one may 
^ ■ resort to an expansion in order to photograph, observe or study it. Actually, expansion of 
the condensate formed on an optical lattice was crucial in the confirmation of Josephson 
oscillation across an array of one-dimensional junctions 4], of Bloch oscillation and 
Landau-Zenner tunneling j^j, of superfluid to Mott insulator quantum phase transition 
and of transition from superfluid to a classical insulator 0E] . Hence it is of interest 
to perform a theoretical study of the expansion of the BEC formed on a joint one- 
dimensional optical-lattice and axially-symmetric harmonic traps upon the removal of 
either trap and compare with available experiments Here we undertake 

this study using the mean-field axially-symmetric Gross-Pitaevskii (GP) equation ^ 
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There have been similar studies of free expansion of an axially-symmetric BEC in the 
absence of an optical-lattice potential [TI] . 

For the present study we consider the expansion of a BEC on a joint one- dimensional 
optical-lattice and axially-symmetric harmonic potentials with the optical-lattice aligned 
along the axial direction of the harmonic potential. First we consider switching off the 
harmonic trap alone. In this case the condensate expands along the radial direction 
maintaining an almost constant transverse size. We also reconsider the free expansion 
of such a BEC formed in the absence of an optical-lattice potential. In both cases the 
time evolution of the e _1 half-width of a Gaussian fit to the density profile in the axial 
and radial directions are obtained and found to be in good agreement with experiments 
by Morsch et al. [TlU lllj. Next we consider switching off the optical- lattice trap alone. 
The periodic one-dimensional optical-lattice potential divides the whole BEC in parallel 
slices. Upon the removal of the optical-lattice trap these slices expand and interfere 
with each other. The result is the formation of an interference pattern of three peaks. 
However, in the presence of the harmonic trap the side peaks execute simple-harmonic 
oscillation with a definite frequency and amplitude. If the interference pattern is allowed 
to evolve in the harmonic trap for some hold time and the harmonic trap is then switched 
off, prominent interference pattern continues for small hold time. With the increase of 
hold time the interference pattern disappears and revives periodically. In this case the 
results of the simulation is in agreement with experiment of Miiller et al. \V2\ . 

In section 2 we present our mean-field model based on the GP equation. In section 3 
we present the numerical results and a comparison with available experiments. Finally, 
the conclusions are given in section 4. 



2. Mean-Field Model 



We base the present study on the numerical solution of the time-dependent GP equation 
[T3*j in the presence of a combined axially-symmetric harmonic and optical-potential 
traps. The time-dependent BEC wave function \I/(r;t) at position r and time t is 
described by the following mean-field nonlinear GP equation ^3] 

-ih^ t -^ + V{r) + gN\^{v-t)\ 2 *(r;t) = 0, (1) 

where m is the mass and N the number of atoms in the condensate, h = h/(2ir) 
with h the Planck's constant, g = ATrh 2 a/m the strength of atomic interaction, a 
the scattering length. In the presence of the combined harmonic and optical traps 
V^(r) = ^mu 2 (p 2 + v 2 z 2 ) + V Q p^ where uo is the angular frequency of the harmonic 
trap in the radial direction p, u z = uoj that in the axial direction z, and V^p^ is the 
optical-lattice trap introduced later. The normalization condition is J dr\^(r; t)\ 2 = 1. 

In the axially-symmetric configuration, the wave function can be written as 
\&(r, t) = ^>(p, z, t). Now transforming to dimensionless variables p = \/2p/l, z = \/2z/l, 
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t = tu, I = y^h/(muj), and (p(p, z; r) = pJ I 3 / \/8ifj(p, z;t), ((TJ) becomes [E] 



.8 d 2 Id d 2 1 




(2) 



where nonlinearity n = Na/l. In terms of the one-dimensional probability 



P(z,t) = 2tt fZ°dp\<p{p,z,T)\ 2 /p 



(3) 



the normalization of the wave function is given by dzP(z,t) = 1. The probability 
P(z, t) is useful in the study of the present problem under the action of the optical- 
lattice potential, specially in the investigation of the formation and evolution of the 
interference pattern after the removal of the trapping potential(s). 

The mean-field GP equation (JTJ) can be derived for a dilute and weakly-interacting 
system. This is based on the condition that the scattering length a be much smaller than 
the average distance between atoms and N » 1. Hence, the dilute gas approximation 
can be written as n\a\ 3 << 1 [13] . where n is the average density of the gas. However, for 
a simple and practical estimate of the weakly interacting condition one should compare 
the interaction energy E- ir ^ with the kinetic energy of the system -E^m- From theoretical 
consideration it follows that 13 j 



The GP equation is valid for small N\a\/l. Systematic previous studies have shown this 
equation to be valid for N\a\/l as large as several hundreds [IH] and the domain of the 
experiments of fOJElE] is well within this limit. 

In the experiments of fDJEJE] a completely asymmetric trap is employed with 
trapping frequencies in three perpendicular directions x, y and z taken in the ratio 
2:l:\/2. An analysis of the experiments using such a fully asymmetric trap is beyond 
present numerical possibilities. However, in the experiment only a section of the 
condensate is observed in one of the x-y, y-z, or z-x planes, e.g., the x-z plane. The 
optical-lattice potential is set in one of the in-plane directions, e.g., the z direction. To 
make the problem numerically tractable, the trapping frequency in the perpendicular y 
direction will be set equal to that in the x direction. This seems to be reasonable as the 
anisotropy in the experiment is not too large. 

We shall consider several aspects of free expansion of a BEC formed on a horizontal 
optical-lattice plus harmonic potentials in the experiments by Morsch et al. [Tn]ITT] and 
Miiller et al. ^2]- in these experiments with repulsive 87 Rb atoms the imaging is done 
in a certain direction and the observed trapping frequencies are ui x = 2n x 25 Hz in 
the radial direction (p) and u z = iov = 2tt x 35.5 Hz in the axial optical-lattice (z) 
direction. The frequency in the unobserved perpendicular direction y is approximated 
by that in the x direction. The actual experimental frequency in the y direction is 
uj y = 2n x 25/ \pl Hz. The optical potential created with the standing- wave laser field is 
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(4) 
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given by V^p^ = VoERsm 2 (nz/d), with Er = h 2 7r 2 /(2md 2 ), Vq the strength, and d the 
lattice spacing. For the mass m = 1.441 x 10~ 25 kg of 87 Rb the harmonic oscillator length 
I = a/ hf (muj) = 2.160 /im and the present dimensionless time unit u; -1 = l/(2ir x 25) 
s = 6.37 ms. In terms of the dimensionless spacing do = \f2djl and the dimensionless 
standing-wave energy parameter Er/(Huj) = n 2 jd\, of is 
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(5) 



Although we employ the dimensionless space units p and z and time unit r in numerical 
calculation, the results are reported in actual units r pm, z pm and t ms and compared 
with the experiments of Morsch et al. JOIE] and Miiller et al. |12j . 



3. Numerical Results 



We solve (J2J) numerically using a split-step time-iteration method with the Crank- 
Nicholson discretization scheme described recently [16] . The time iteration is started 
with the harmonic oscillator solution of (j2J with n = 0: ip(p,z) = [u / '(87T 3 )] 1 / 4 
p e -(p 2 +^^ 2 )/4 First the nonlinearity n is slowly increased by equal amounts in lOOOOn 
steps of time iteration until the desired value of n is attained. Then the optical-lattice 
strength Vq is increased in 10 000 small steps and the time iteration continued until the 
final value of Vq is attained. Then, without changing any parameter, the solution so 
obtained is iterated 50 000 times so that a stable solution is obtained independent of 
the initial input and time and space steps. 

The BEC on the optical lattice for a specific nonlinearity and the interference 
pattern upon the free expansion of such a BEC have been recently studied using the 
numerical solution of the GP equation while both the optical-lattice and harmonic 
potentials are switched off in one, two and three dimensions 0HZI- Here we study 
the expansion of the BEC when either the harmonic potential or the optical-lattice 
potential is switched off. The periodic optical-lattice potential cuts the condensate into 
several pieces in different sites. As the present calculation is performed with the full 
wave function without approximation, phase coherence among different optical-lattice 
sites is automatically guaranteed in the initial state. As a result when the condensate 
is released from the periodic optical-lattice trap, a matter-wave interference pattern is 
formed in a few milliseconds. The atom cloud released from one lattice site expand, 
and overlap and interfere with atom clouds from neighboring sites to form the robust 
interference pattern, consisting of a central peak and two smaller symmetrically spaced 
peaks moving in opposite directions. If the harmonic potential is kept on while switching 
off the optical-lattice potential the two side peaks execute simple harmonic oscillation 
with the axial frequency uj z . This simple harmonic motion can be explained considering 
the dynamical evolution of the interference peaks in the harmonic potential [T2"lll8j . 
Since the lattice transfers momentum to the condensate in units of 2p& = h/d, the 
recoil velocity is given by 2vr = 2pn/m = h/(md) and the two side peaks move with 
velocities ±h/(md) at time t — immediately after release. If the harmonic trap is 
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Figure 1. View of the condensate in the radial direction after an expansion time 
^exp = 0, 5 ms, 10 ms, 15 ms, 20 ms and 25 ms in the optical-lattice trap alone 
obtained from the contour plot of the wave function \ip(p, z, £exp)| for Vq — 25, 
n = 53, uj = 2tt x 25 Hz and lo z = 2tt x 35.5 Hz. To obtain the full section we 
used \i/j(-p,z,t exp )\ = \ip(p,z,t e xp)\- 

maintained the two peaks will oscillate with an amplitude h/{mduj z ). Consequently, the 
motion of the side peaks is given by [T2"] 

z(t) =±— J— sin(u z t). (6) 

The experiments of Morsch et al. fTUlITT] and Miiller et al. [12] were performed with 
about (1 — 3) x 10 4 Rb atoms with the optical lattice spacing d = 1.56 jum. In our 
simulation we use N = 20000 atoms, to = 2ir x 25 Hz and uj z = 2n x 35.5 Hz. With the 
scattering length a = 5.77 nm, this corresponds to a nonlinearity of n = Na/l « 53.4. 
First we consider the results of free expansion reported in figures 1 and 2 of [TU] and 
figure 1 of JT]. These results are related to the expansion of the condensate when only 
the harmonic potential is switched off. Later we consider the results reported in ^2] 
about the simple harmonic oscillation of the interference peaks of the condensate when 
only the optical-lattice potential is switched off as well as the disappearance and revival 
of interference pattern after a time evolution in the harmonic trap. 

To investigate the expansion of the condensate when the harmonic trap alone is 
switched off we use optical-lattice strength Vq = 25. Once the condensate is prepared 
on the joint optical-lattice and harmonic potentials the harmonic potential is switched 
off and the expansion studied under the influence of the optical trap alone. The view 
of the condensate at different expansion times texp = 0, 5 ms, 10 ms, 15 ms, 20 ms 
and 25 ms are shown in figure 1. The view at (20 — 25) ms is similar to figure 1 (b) 
of JU] after 23 ms. The strong optical lattice divides the condensate into smaller pieces 
trapped in different lattice sites. With time one piece of BEC in a particular site expand 
perpendicular to the lattice without ever moving to a neighboring lattice. Consequently, 
the BEC expands only in the radial direction and maintains a constant size in the axial 
direction. 




Figure 2. Condensate dimensions (p) (radial) and (z) (axial) versus the expansion 
time iexp in the presence of a lattice of strength (a) Vq = 25 and (b) Vo = 0. The 
nonlinearity and frequencies are the same as in figure 1. Open square: (p) of the 
experiment of |TJ3]; star: (z) of the experiment of JUj; full line: (p) of the present 
simulation; dotted line: (z) of the present simulation. 
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Figure 3. The ratio (z)/(p) versus the expansion time texp in the presence of lattices 
of strength Vo — 20 and Vb = 0. The nonlinearity and frequencies are the same as 
in figure 1. Open square: experiment of ^ for Vq = 0; star: experiment of for 
Vo = 20; dotted line: present simulation for Vo = 0; full line: present simulation for 
Vo = 20. 

Next we calculate (z) and (p) which are the e _1 half-widths of Gaussian fits to 
the density profile in the lattice (axial) and radial directions, respectively, for different 
expansion times. The results are plotted in figure 2 (a) for Vo = 25 and compared with 
the experiment of JU] ■ in figure 2 (b) we plot the same results of expansion as in figure 
2 (a) for Vo = and compare with the experiment ^0] • Considering the 5 fim resolution 
of the imaging system ^T] the agreement between the present numerical simulation and 
experiment is good. 

Next we perform the same expansions of figure 2 with optical lattices of strength 
Vq = 20 and and plot the ratio (z) / (p) in these cases for different expansion times and 
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Figure 4. Contour density plot of one-dimensional probability P(z, t) of © versus z 
and £ for the initial BEC of figure 1 after the removal of the optical trap at t = when 
an interference pattern of a central peak and two side peaks is formed. The sinusoidal 
pattern of oscillation appears in the contour density plot and is compared with the 
simple theoretical prediction ©, e.g., z(t) — ±13.3sin(27r x 35.5i) /im, shown by the 
dark solid line. 



compare with the experimental results of llj. The agreement is good. In both figures 
2 and 3 the agreement is better for large expansion times. At large expansion times 
the size is larger and hence possibly has less experimental error thus leading to a better 
agreement with experiment. 

Miiller et al. [T2| performed a different type of expansion experiment on a BEC 
formed on a joint optical-lattice and harmonic traps upon switching off only the optical- 
lattice trap. The condensate pieces from different sites expand, overlap and interfere to 
form a pattern with three peaks which oscillate according to ©• In figure 4 we show 
the contour density plot of P(z, t) of (jHJ) as functions of z and t. This clearly shows 
the simple harmonic motion of the two side peaks which is compared with the simple 
theoretical prediction (JBj) . The agreement between the two is good. 

Finally, Miiller et al. JE] observed in a remarkable experiment the evolution of the 
interference patterns (formed upon switching off the optical-lattice trap) in the harmonic 
trap alone for different times. Then they removed the harmonic trap and observed the 
free expansion of the BEC. For small hold times T in the magnetic trap the interference 
pattern is clearly visible. They found that with the increase of hold time the prominent 
interference pattern first becomes blurred and then disappears completely. Then it 
reappears prominently after a complete half cycle of oscillation according to ©• In 
this case the period of oscillation is 28.2 ms. The numerical simulation at t = 0.1 
ms in figure 5 (a) shows the clear trails corresponding to the interference pattern in 
agreement with [I2j- The complete free expansion for a hold time T = leads to a 
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Figure 5. Contour density plot of one-dimensional probability P(z, t) of J3J) versus 
z and t for the initial BEC of figure 1 after the removal of the optical trap at t = 
when an interference pattern of a central peak and two side peaks is formed. The BEC 
is then kept in the harmonic trap for hold times (a) T = 0.1 ms, (b) T = 6 ms, (c) 
T = 10 ms and (d) T — 18 ms before finally removing it. The BEC is then allowed 
to expand for 15 ms in order to see the interference peaks. The interference pattern 
of three peaks represented by three trails appears for small hold times (a), becomes 
blurred (b) and disappears (c) with the increase of hold time and finally reappears (d) 
when the hold time is larger than the half period of oscillation of the peaks given by 

©■ 
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plot indistinguishable from figure 5 (a). In this case the velocity of the side peaks is 
±h/(md) = 2.96 /xm/ms. After 15 ms of free expansion each side peak moves 44.4 /jm 
consistent with the simulation of figure 5 (a). In figure 5(b) it is blurred at t = 6 ms 
in agreement with ^2]- The interference pattern is completely destroyed at t = 10 ms 
and again revived after a half cycle of oscillation at t — 18 ms in figures 5 (c) and (d), 
respectively, in agreement with experiment. 

The explanation for the phenomenon observed in figures 5 is as follows. As the 
system oscillates in the harmonic potential upon the removal of the optical-lattice 
potential, after every half cycle of oscillation the system passes through the same initial 
state when the optical-lattice potential is removed fZj. This initial state is a coherent 
state on the optical lattice and after a free expansion leads to the interference pattern. 
However, if the system is allowed to evolve in the harmonic potential for some time after 
removing the optical-lattice potential and then subject to a free expansion a completely 
different scenario emerges. In general, the three expanding bright spots would now be 
moving freely in "wrong" directions so that they may interact and collide with each 
other upon removal of the harmonic trap and lose their identity resulting in a loss of 
the observed interference pattern in figures 5. This defocussing effect in the interference 
pattern would be absent for hold times equal to a multiple of half period of oscillation 
and would be minimal for hold times near these values consistent with figures 5. There 
has been a recent attempt to study a similar phenomenon in a different context using a 
one-dimensional mean- field model |19j . 

4. Conclusion 

In conclusion, using the explicit numerical solution of the axially-symmetric GP equation 
we have studied the expansion of a BEC formed in a joint harmonic and optical-lattice 
traps when either trap is removed. When the harmonic trap is removed the piece of 
condensate in a particular site expands radially without moving to the next site. As a 
result the whole condensate expands radially maintaining a fixed axial size. The shape 
of the condensate and its radial and axial sizes during expansion are found to evolve in 
agreement with the experiments of [TTHlll| . When the optical trap is removed the pieces 
of the BEC at different sites expand and interfere to form a distinct pattern with three 
peaks which keep on oscillating according to (jBJ) in agreement with present simulation 
and experiment of [T2"] . If the interference pattern is allowed to evolve in the harmonic 
trap for some hold time and the harmonic trap removed, prominent interference pattern 
appears for small hold time, which become blurred and finally, disappear with the 
increase of hold time in agreement with the experiment of . Finally, for hold times 
larger than half period of oscillation the interference pattern reappears as in experiment. 
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